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Abstract. We investigate properties and the distribution of light nuclei (A < 4) in symmetric nuclear 
matter of finite temperature within a microscopic framework. For this purpose we have solved few-body 
Alt-Grassberger-Sandhas type equations for quasi-nucleons that include self-energy corrections and Pauli 
blocking in a systematic way. In a statistical model we find a significant influence in the composition of 
nuclear matter if medium effects are included in the microscopic calculation of nuclei. If multiplicities 
are frozen out at a certain time (or volume), we expect significant consequences for the formation of 
light fragments in a heavy ion collision. As a consequence of the systematic inclusion of medium effects the 
ordering of multiplicities becomes opposite to the law of mass action of ideal components. This is necessary 
to explain the large abundance of a-particles in a heavy ion collision that are otherwise largely suppressed 
in an ideal equilibrium scenario. 

PACS. 25.70.-z Low and intermediate energy heavy-ion reactions - 25.70.Pq Multifragment emission and 
correlations - 21.65.+f Nuclear matter - 21.45.+V Few-body systems 
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1 Introduction 

Heavy ion collisions provide a tool to investigate the phase 
structure of nuclear matter. Depending on the energies, 
the region of temperature and density explored might be 
rather large. The information about the composition of 
nuclear matter is contained in the equation of state. At col- 
lision energies per nucleon well below one GeV the equa- 
tion of state is described by purely hadronic degrees of 
freedom. It is a basic ingredient in microscopic simulations 
of the heavy ion collision, such as the Boltzmann-Uehling- 
Uhlenbeck (BUU) [rU^lTflliHB"] or the quantum molecular 
dynamics (QMD) |^]|7||B) simulations. The challenge is to 
extract information about the different stages of the evo- 
lution of the heavy ion collision. This information could 
be provided by fragments produced in the different stages 
of the collision as has been done recently for the case 
129 Xe+ nat Sn by the INDRA collaboration |3ITcT] . 

An early analysis of multi-fragmentation in a heavy ion 
collision of 36 Ar+ 58 Ni at several energies below 100^4 MeV 
has been given in Ref. The authors study a class of 
evaporation events at central collisions |12|. During these 
events more than 90% of the charged particles were de- 
tected and isotopically identified. Within a thermal (and 
chemical) equilibrium scenario |KSj of ideal gas compo- 
nents (including states up to excited 9 B), supplemented 
by finite volume effects jH] and a model of side-feeding, 
they found a remarkable agreement with the experimen- 
tal data Temperature has been varied between 10 and 



25 MeV and the freeze-out volume fixed to 1 /3 of normal 
nuclear matter density. 

A more elaborated statistical analysis has been done 
for the recent INDRA experiment 129 Xe+ nat Sn. The mea- 
sured multiplicities of the central collisions by the INDRA 
collaboration show a large fraction of a-particles jni|. In 
contrast, a naive model of a gas of ideal components, 
would give a much smaller number of a-particles, depend- 
ing on the freeze-out density. The INDRA collaboration 
provides a detailed comparison of their data within a sta- 
tistical multi-fragmentation model (SMM) . This model 
goes beyond a simple picture of an ideal gas and describes 
multiplicities and some other aspects of the heavy ion col- 
lision in question |l()j . 

On the other hand, from a microscopic analysis of clus- 
ter formation, it is known that nuclei dissociate already at 
rather moderate densities and temperatures, see, e.g. ^3 
IT7] and references therein. Some details will also be given 
in this paper. The dissociation (Mott effect) is taken into 
account, e.g., in modern BUU simulations of heavy ion 
collisions and is necessary to reproduce the experimental 
data, see Ref. It is an effect related to the Pauli 

blocking induced by the surrounding medium and goes 
beyond the picture of a simple ideal gas of nuclei. 

In a recent analysis of the central collision Xe+Sn at 
50A MeV that has been measured by the INDRA collabo- 
ration [2] , we found that the BUU simulation gives a pro- 
ton to deuteron ratio that is close to the one expected by 
the equilibrium distribution; to be more precise, for times 
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t > 50 fm/c during the evolution of the system. This how- 
ever holds only, if the equilibrium distribution includes the 
above mentioned dissociation of the deuteron [olll8| . 

Therefore, we address the question to what extend the 
dissociation of nuclei affects the equilibrium distribution 
of nuclear matter. To do so, we investigate a system of 
light nuclei at finite density and temperature up to the 
a-particle. 

Hence, we focus on a new aspect in the distribution of 
light nuclei. Since reasonable generalizations of the Feynman- 
Galitskii or Bethe-Goldstone equations for more than two 
particles become available, the properties of multi-particle 
correlations in a medium can now be addressed micro- 
scopically. To demonstrate the effect that is related to 
in-medium properties of the light clusters in question, in 
particular the a-particle, we explore an ab initio equilib- 
rium quantum statistical description of a many-particle 
system based on the well established and successful (equi- 
librium) Green function method . To include a proper 
description of clusters we implement an equal time con- 
straint on the Green functions. This allows for a cluster 
expansion of the Green functions as shown, for example 
in |20|. In an uncorrelated medium of quasi-particles the 
equal time constraint systematically leads to Dyson equa- 
tions for clusters with a fixed number of particles. They 
include the self energy corrections and the Pauli block- 
ing and are rearranged as resolved equations to use the 
Alt-Grassberger-Sandhas (AGS) formalism to solve the 
respective few-body equations which has been done for 
the three- and four-nucleon system in |21U16II17| . Similar 
equations to treat the in-medium three-body system have 
been proposed previously in Refs. [23124] . As an interac- 
tion we use a nucleon nucleon potential that reasonably 
reproduces the nucleon nucleon phase shifts and the bind- 
ing energies of the light nuclei in question. 

In Section[5]we introduce the consequences of the above- 
mentioned cluster expansion method in the equation of 
state. This will be done along the lines of The mi- 
croscopic AGS-type equations to treat multi-particle clus- 
ters in medium will be explained in Section |3| We use a 
nucleon-nucleon potential that reproduces reasonably well 
the nucleon-nucleon scattering data and the binding ener- 
gies of light nuclei considered. In Section^we present our 
results. In particular, we calculate the equilibrium compo- 
sition of nuclear matter for conditions comparable to the 
heavy ion collision investigated by the INDRA collabora- 
tion. We summarize our conclusions in Sectional 



2 Statistical model 

To generalize the equation of state for a Fermi system 
that includes correlations the nuclear matter density n = 
n(/x, T) as a function of the chemical potential fi and tem- 
perature T can be rearranged in an uncorrelated part nf ree 
and a correlated one n corr 25,.2j, 

71 ^frcc ^corr- (I ) 



To abbreviate notation let 1 denote the quantum numbers 
of particle 1. The Fermi function is given by 

f(l) ee f{E x ) = {exp^ - ,i)] + l}- 1 (2) 

where E\ denotes the one-particle energy and (3 the in- 
verse temperature. Presently we describe symmetric nu- 
clear matter and hence 

n free = 4^/(1). (3) 
i 

For a system of nucleons of mass rajy the energy is E\ = 
k 2 /2mpf. Hartree-Fock approximation introduces the no- 
tion of quasi-particles and quasi-particle energies 

E x ->£i = k 2 /2m N + £(k), (4) 

where 

£(l) = $> 2 (12,12)/(2). (5) 

2 

The tilde means proper anti-symmetrization (i.e. includ- 
ing the Fock term). 

As explained in Ref. ^Hj m some detail, the correlated 
density can be composed into different cluster contribu- 
tions. 

^corr = 2n 2 + 3n 3 + 4n 4 + . . . , n A — n A + n A (6) 

where n A denotes the A-particle correlated density pre- 
sented as bound n\ or scattering states in chemical 
equilibrium. The full expression for the two-particle cor- 
related densities n 2 has been given in Ref. |25| . 

To evaluate correlated densities we presently focus on 
the bound state contributions. This is justified in view 
of the rather low densities of the final stage of heavy ion 
collisions. Note, however, that including contributions of 
scattering states requires major theoretical effort to solve 
the respective scatting few-body equations derived in the 
next section. The distribution functions for the A-body 
cluster of fermions is given by 

f A (p) = {exp\f3(E A -B A -ti A )]+e}- 1 (7) 

where p is the cm. momentum of the cluster, Ea{p) is 
the continuum energy, B A > the binding energy of the 
cluster, e = +!(—!) for fermions (bosons), and ha the 
respective chemical potential. In equilibrium considered 
here \i A = A/j,. The density for the nucleus of mass number 
A is given by 

n A ^T) = (2S + l)(2I+l)J2fA(p) (8) 

p 

where S denotes the spin and / the isospin of the nucleus. 



3 In-medium few-body equations 

The basis of the equations given in this section is the (equi- 
librium) Green function approach to describe quantum 
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statistical systems ^Hj- The Green functions for a given 
number of particles are evaluated at equal imaginary times 
assuming an environment of independent quasi-particles, 
see e.g. [20]. As a consequence, for a given number of 
particles, a Dyson-type equation can be derived that is 
only driven by the dynamics of the smaller cluster (clus- 
ter mean field approximation) that breaks the Green func- 
tions hierarchy. Utilizing resolvents the Dyson-type equa- 
tion for a particular cluster can be rewritten as AGS-type 
equation 26,27,28, 29 with an effective Hamiltonian. This 
has been shown previously for the nucleon deuteron reac- 
tion [21] the three-nucleon bound state ^Hj an d the a- 
particle ^Zj ■ In this section we briefly repeat the relevant 
formulas to introduce our notation. 

Utilizing the Dyson equation for clusters it is possible 
to introduce resolvents to describe the dynamics of the 
system. Defining H = Y^i=i s i-> w ^ n the quasi-particle 
self energy the n-quasi-particle cluster resolvent Go is 

GoO) = {z - Hq)- 1 N = Rq(z) N. (9) 

Here Go, Hq, and N are matrices in n particle space and 
z denotes the analytic continuation of the Matsubara fre- 
quency JU] . The Pauli-blocking factors for n-particles are 

N = /(l)/(2) . . . /» + e/(l)/(2) . . . f(n), (10) 

with / ee 1 - /. Note: NR = R N. Defining the effec- 
tive potential V = Epairs a N%Vg the full, G(z), and the 
channel, G a (z), resolvents are 

G(z) = (z - H - V^N ee R(z)N, (11) 
G a (z) = (z-H ~ N%V?)- l N ee R a {z)N. (12) 

Note that U 1 * ^ V and R(z)N ^ NR(z). For the scatter- 
ing problem it is convenient to define the in-medium AGS 
operator Up a (z) |2J 

R(z) = Sp a R (z) + R p {z)U 0a {z)R a {z) (13) 

that after some algebra leads to the in-medium AGS equa- 
tion 

Up a (z) = E^Roiz)- 1 +Y,S f 3- ( N^(z)Ro(z)U ia (z), 

7 

where 8p a == 1 — 5p a . The square of this AGS-operator 
is directly linked to the differential cross section for the 
scattering process a — > /3, for all Fermi functions f(i) — > 0. 
Hence the isolated three-body system is recovered. The 
driving kernel consists of the two-body t-matrix derived in 
the same formalism, however given earlier and known as 
Feynman-Galitskii (finite T) or Bethe-Goldstone (T = 0) 
equations |19II22| 

A numerical solution of the three-body break-up reaction 
relevant for the chemical distribution in a heavy ion colli- 
sion using a coupled Yamaguchi potential has been given 
in Ref. [2IJ. 



For the bound state problem it is convenient to intro- 
duce form factors 

i^> = EWWi^3>- (15) 

7 

Since the potential is non symmetric, the right and left 
eigenvectors are different, although the bound state ener- 
gies are the same. The eigenvectors are explicitly needed 
in our solution of the four-body system. The respective 
homogeneous AGS equations are given by 

\F a ) =J2^N^Ti(B 3 )Ro(B 3 )\F ), 



\F a ) =J2^Ti(B 3 )NiMB 3 )\F ), (16) 



We now turn to the four-body problem in matter. In 
addition to having different channels as for the three body 
system now the channels appear in different partitions 
that makes the four-body problem even more involved. 
The partitions of the four-body clusters are denoted by 
p,r,a,..., e.g., p = (123)(4),(234)(1),... for 3 + 1-type 
partitions, or p = (12)(43), (23)(41), ... for 2 + 2-type par- 
titions. The two-body sub-channels are denoted by pair 
indices a, 0,-y, . . ., e.g. pairs (12), (24), ...The two- and 
three-body ^-matrices have to be defined with respect to 
the partitions that leads to additional indices. A conve- 
nient way to solve the four-body in-medium homogeneous 
AGS equation is by introducing form factors 

i^) = E^E^ i? o 1 (s 4 )i^ 4 ), (17) 

T a 

where 5J ja — Sp a , if /3, a C r and 5p a = otherwise 
and \ipBi) is the a-particle in-medium wave function. The 
homogeneous equations then read |17| 

(18) 

where a C er, 7 C r. A numerical solution of this equa- 
tion is rather complex. In order to reduce computational 
time, needed in particular to handle the dependence on 
the medium, we introduce an energy dependent pole ex- 
pansion (EDPE) that has been proven useful in many ap- 
plications involving the a-particle and is accurate enough 
for the present purpose |3U| . However, we have to general- 
ize the original version of the EDPE because of different 
right and left eigenvectors appearing for the three-body 
subsystem and given in Eq. l|16f) (for details see |17|). 
In the two-body sub-system the EDPE reads 

T li Z ) — E \ r in{z))t ln (z)(r in (z)\ 
n 

— ^ ^ \9'yn)t-fn( z ){gyn\ 
n 

= Y, N 2\9m)hn{z){g in \. (19) 
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where we have chosen a Yamaguchi ansatz for the form 
factors for simplicity. The last line shows the explicit de- 
pendence of the Pauli blocking factors. Inserting this ansatz 
into the Feynman-Galitskii equation determines the prop- 
agator t in (z). In the three-body sub-system the EDPE ex- 
pansion reads 

{g(3rn{z)\Rn{z)U T M (z)R Q {z)\g in (z)) ~ 
with the three-body EDPE functions 

= (9*n\Mz)\hrn)tp m (B S )\f;^), (21) 

that we get from solving the following Sturmian equations 
Vtjf^) = Y.(9 a n\R O {B^\~g 0m )t Pm (B^)\r;^) (22) 

0m 

mjr^) = J2(9 a n\Ro(B3)\g Pm )t Pm (B 3 )\r^) (23) 

0m 

Inserting the EDPE into the homogeneous AGS equa- 
tions allows us to redefine the form factors that are now 
operators depending on the coordinates of the 2 + 2 or 

3 + 1 system, i.e. 

IF?*) = E^,,(^4)IW( S 4)(^m(54)|Eo(54)|J^) 

(3m 

and therefore the final homogeneous equation 

rt uk 771 

xt.M) \r?), (25) 

is an effective one-body equation with an effective poten- 
tial V and an effective resolvent Qq defined as 

VP T( W =J2^r(r^(z)m n (z)\f^(z)}, (26) 

(27) 

4 Results 

4.1 Cluster properties 

The binding energies of the few-nucleon systems depend 
on the chemical potential /i or equivalently the density 
n(/i,T), the temperature T, and cm. momentum P c . m .. 
For the two-, three-, and four-nucleon systems the bind- 
ing energies are shown in Fig. for T = 10 MeV and 
-Pc.m. = |T7]. The line Ba — reflects the respective 
continuum threshold. We mention here that the medium 
dependence of the binding energies is rather similar for dif- 
ferent two-body potentials, although their results for the 
isolated system may be very different for the few-nucleon 
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Fig. 1. Density dependence of the binding energy per nucleon 
of deuterons (dash-dot), Triton (dashed), and a-particle with 
Malfliet Tjon potential (solid), Yamaguchi potential (long- 
dashed) at Pc.m. = 0, as given in Ref. |17|. 
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Fig. 2. Momenta per nucleon as a function of the dissocia- 
tion density for deuteron (dash-dot), triton (dash), a-particle 
(solid). Bound states exist only above the respective lines. 



systems considered. This has been mentioned earlier, but 
it can also be seen from the two lines representing dif- 
ferent potentials for the a-particle, i.e. Yamaguchi |3J 
(long dashed) and Malfliet-Tjon potential [S3 (solid) after 
renormalizing the binding energies to the same value of the 
MTI-III potential. For 3 He that is not shown in Fig.Q]thc 
dissociation density is slightly smaller due to the Coulomb 
force that has been evaluated perturbatively. However, for 
asymmetric nuclear matter, e.g. N p /N n ~ 0.72 (for the 
129 Xe+ 119 Sn reaction) this effect is compensated |33| . 

Because of the medium the pole of the bound state 
moves to the continuum threshold as seen in Fig. ^ The 
bound state vanishes because of the quasi-particle na- 
ture of the cluster. However, investigating the zeros of the 
two-body Jost function we found earlier that the quasi 
deuteron "survives" as a virtual bound state with differ- 
ent energies, depending on the densities above the dis- 
sociation line This is a similar state as the virtual 
1 S'o nucleon nucleon state at 70 keV. This means that 
the quasi deuteron retains its quasi-particle nature of be- 
ing an infinitely long living state and does not become a 
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0,15 



Fig. 3. Density dependence of the binding energy of the a- 
particle for finite P c . m . as indicated in the legend. 



For a given temperature, here we chose T = 10 MeV, 
the momentum of dissociation Pdis is defined by the con- 
dition 

B(n dis ,T,P dis ) = 0, (28) 

i.e., the density n^is and the momentum Pdis where bind- 
ing of the nucleons is lost. For a system of atoms and ions 
this scenario can be related to a transition of an isolator to 
a conductor, since electrons can move away from the ion 
because of the dissociation property. The respective dis- 
sociation lines for deuterons, tritons/ 3 He, and a-particles 
are shown in Fig. [5] The momentum is normalized to the 
number of nucleons in the cluster that is identical to the 
velocity of the cluster on the dissociation line. 

4.2 Cluster distribution 



resonance. Only going beyond the quasi-particle picture 
the deuteron spectral function will aquire an imaginary 
part and hence the deuteron becomes a state of finite life 
time in the medium. This is due to break-up processes in 
nuclear matter, the simplest one being the three-nucleon 
reaction 12 1| . Such an investigation for three- and four- 
body system is also technically involved and still needs to 
be done. 

Also a possible appearance of Efimov states related to 
B — > of the sub-system needs further investigation |35j . 
Since the Efimov states are 'excited' states, e.g. for the 
three-body system they are close to the 2 + 1 threshold, 
their blocking may be smaller since the wave functions 
contain higher momentum components, hence the slope 
of their dependence on the density is flatter. Also, as seen 
from Fig. Q the slope of the binding energies as a func- 
tion of densities for the larger clusters is steeper. On the 
other hand the sub-system is not at rest in the larger clus- 
ter, hence the binding energy changes, as we will show 
in the next paragraph, and therefore a careful analysis 
is needed that would go beyond the present scope of the 
paper. Hence so far no conclusion can be drawn for the 
appearance of Efimov states, but it is an important issue 
since Efimov states might effect the equation of state for 
clustering Fermi systems. 

For a finite cm. momentum relative to the medium (at 
rest) the influence of the medium is weaker, as less com- 
ponents of the wave functions are blocked by the Fermi 
sea. For deuteron j^H] and triton ^1 this has been given 
in earlier references, see also references therein. In Fig. |3 
we give the results for the a-particle. Note again that the 
medium effects do not change the elementary property of 
an a-particle, however, after introducing effective degrees 
of freedom, the a-particle and any other cluster considered 
here consists of quasi-nucleons and not elementary nucle- 
ons. Besides the change of nucleon self energy also the 
binding energy of the cluster is changed and hence the 
clusters can be viewed as quasi-deuterons, quasi-tritons, 
quasi-a's, etc., i.e. clusters with the respective self-energy 
corrections. 



We now consider the composition of nuclear matter in 
equilibrium at a temperature of T = 10 MeV. To this end 
we assume that nuclear matter is composed of nucleons, 
deuterons, tritons, 3 He, and a-particles. Larger clusters 
are presently not considered. We investigate three scenar- 
ios: 

A. gas of nucleons and nuclei with all properties of the 
isolated systems retained (no change due to medium), 

B. a gas of quasi-nucleons that contain self energy cor- 
rections, the clusters consist of quasi-nucleons but the 
interaction is without Pauli blocking, hence the bind- 
ing energy will not change, and 

C. a gas of quasi-nucleons with Pauli blocking in the inter- 
action and therefore the clusters are treated as quasi- 
nuclei that include the self energy corrections on the 
cluster level and dissociation. 

In the ideal situation (case A), the components retain their 
properties, i.e. all particle masses stay the same. The com- 
position of the system is driven by the law of mass action, 
i.e. the equilibrium distribution functions of nuclei (con- 
sisting of A nucleons) with mass m A = Attin — B A are 
given by 

f A (p) = {exp[/3(p 2 /2m A - B A - ^ A )} + e}^ . (29) 

The composition of the system, i.e., the number of par- 
ticles as a function of the total density, for T = 10 MeV 
is given in Fig. The density is accumulated by the 
more massive clusters the larger the total densities gets. 
The freeze-out distribution could be read of at n to t = 
0.085... 0.034 fin" 3 Q31- Fig. 0)3 refers to the result of 
the quasi-particle approximation (case B) for nucleons, 
instead of using ideal nucleons. In this case the medium 
effects are taken into account in Hartree-Fock approxima- 
tion on the single particle level. This results in different 
self energies for the nucleon and in turn the mass of nuclei 
changes accordingly, see Eq. 

e(k) ~ k 2 /2m cS + E(0) (30) 

The right hand side of Eq. I|30|) , known as effective mass 
approximation, is valid for the rather low momenta and 
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The situation changes drastically for case C, if the 
change of the binding energy as discussed above is taken 
into account, 



Fig. 4. Numbers of nucleon (long dashed), deuteron (dotted), 
triton ( 3 He) (dashed), 4 He (solid) as a function of the total 
density n to t at T = 10 MeV. Total number of nucleons 250. A) 
ideal system B) quasi particle system, B) full calculation with 
dissociation. 



densities considered, hence m eS = rriN^fi, T) approximately 
independent of the momentum k. The distribution func- 
tions change to 

f A (p) = {exp[/V/2™°/ -B A - (if)] + e}' 1 (31) 



and now 



,cff 



= Am e 



B A . In chemical equilibrium 



(i A — Afi cS , where (l cS = (x — S(0). The results are rather 
close to the ideal gas case, because the change of the self 
energy of the cluster due to the binding energy of the 
bound nuclei is not taken into account. Differences appear 
at larger densities. 



— DCff 



B A ^B A {p,T,n) = B 



(32) 



This, however, needs a solution of few-body in- medium 
equations as given in the previous section. The effects 
induced by this change in binding energy is shown in 
Fig-HP- The equation for the density of the cluster changes, 
because the bound state exists only above the momentum 
of dissociation as shown in Fig. [3 The definition of the 
density changes accordingly, see Eq. JSJ 

r&(/i,T) = (2£ + l)(2J + l) £ f A (p). (33) 

P>Pdia 

Also the distribution function is different from the previ- 
ous definition, since now the change of the binding energy 
has to be taken into account. 



f A (p) = {exp[/3(p 2 /2mf - Bf - (if)] + e}' 



(34) 



where now mf = Arojv(/J, T) — Bf and Bf is given in 
Fig. |3 for the a-particle, for the three-body case in 
and for the deuteron in [3S], and Refs. therein. 

The change between the ideal (or quasi-particle) pic- 
ture and the full calculation that includes self energy cor- 
rections and dissociation of the clusters appears quite de- 
cisive. Whereas for rigid nuclei the number of heavy parti- 
cles is much higher than the number of light particles, this 
is different, if the dissociation is taken into account. The 
geometrical interpretation of dissociation is because less 
low momentum components (large distance components) 
are available for the formation of a bound state. If the 
momentum is higher, the nucleus moves out of the Fermi 
sphere of the surrounding matter and the particle becomes 
more stable. The fact that no bound states are possible 
doesn't mean that there are no correlations. Previously 
we found by analyzing the Jost function of the deuteron 
that as the deuteron moves towards lower binding ener- 
gies and eventually crosses the continuum line (B 2 = 0) it 
exists as a virtual bound state on the unphysical energy 
sheet |34| . In its turn this means that particular correla- 
tions in the continuum can form deuterons, if the density 
becomes low enough (for a given temperature). A simi- 
lar study for three- and four-body states still needs to be 
done. On the other hand scattering states are infinitely 
extended, in contrast to bound states. So neglecting those 
few-body correlations related to scattering states may not 
be the worst approximation to start with. 

To give an example for possible consequences of this 
finding, we investigate the multiplicities of light fragments 
in a heavy ion collision. We investigate conditions close to 
the INDRA experiment Xe on Sn at 50 A MeV jHj studied 
earlier within the context of a BUU simulation We 
focus on the final stage of the collision and for simplic- 
ity assume a homogeneous temperature of T = 10 MeV 
that might slightly be too high for a quantitative compari- 
son, but is still reasonable. Also, we use symmetric nuclear 
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Fig. 5. Evolution of the density in the rrz-plane as provided by a BUU simulation |18| of a central collision of Xe on Sn. 



matter that is not achieved in the experiment mentioned. 
However, the basic effects that we focus on in this paper 
are not changed: We found previously that asymmetric 
nuclear matter (using the asymmetry induced by the ex- 
periment Xe on Sn) has very little effect on the dissocia- 
tion of 3 He and triton. Quantitatively, this effect is in the 
same order of magnitude as the Coloumb correction pTH] , 
The BUU simulation of the central collision of Xe on 
Sn at 50A MeV provides a realistic nuclear density dis- 
tribution for the INDRA experiment. A cut through the 
xz-plane is shown in Fig.l5l[T%]. To simplify and model this 
density distribution we assume a homogeneous spherical 
distribution of radius R that approximately matches the 



size of the simulation at about 40 fm/c (R = 7.5 fm) and 
about 140 fm/c (R = 20 fm) after the collision. The ra- 
dial change is assumed to be linear. The resulting change 
of the local density with respect to time for this simple 
expanding fire ball is shown in Fig. 



The multiplicities of nucleons, deuterons, 3 He/ 3 H, and 
a-particles for the densities evolution of Fig. and a tem- 
perature of T — 10 MeV for a total number of nucleons 
of 250 is given in Fig. [7| This is the main result of the 
present calculation. 
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Fig. 6. Evolution of total density in a uniform model. Param- 
eters are chosen as explained in the text. 
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Fig. 7. Total numbers of nucleons (long-dashed x 10), deuteron 
(dots), triton/ 3 He (dashed), and a particles (solid) as a func- 
tion of collision time at T = 10 MeV. 



4.3 Discussion 

At around t = 100 fm/c the number of a-particles is much 
larger than the number of other clusters. This time, which 
(approximately) corresponds to about 3-4 times the initial 
volume, can be considered as freeze-out time (related to 
the freeze-out volume fTO]). Hence the multiplicities at this 
time should be significantly correlated with the experi- 
mentally observed ones. Indeed, the total multiplicities of 
a-particles in the above mentioned INDRA experiments 
Xe+Sn is larger than those of the lighter clusters ^Oj • 

In contrast to the full calculation, the other equilib- 
rium scenarios discussed above have the opposite order- 
ing of multiplicities which can be concluded from Fig. rj] 
for ntot — 0.04fm -3 and does not reflect the experimental 
finding for the large excess of a-particles. We argue that 
the enhancement of a-particles is also related to the fact 
that the a-particle is more stable in low density nuclear 
matter than the other light clusters. Our findings of the or- 
dering of multiplicities for the light fragments obtained via 
a microscopic approach and including dissociation might 
provide a natural explanation to the excess of a-particles. 
However, for a throughout comparison with experimental 
data several other aspects have to be taken into account 



as mentioned in jlOlllr?) . However, some of them need a 
major effort while going beyond the model of an ideal gas 
of components. 

Concerning larger (light) clusters than the ones con- 
sidered so far, note that they are weaker bound than the 
a-particle. Therefore they should be less stable in medium. 
Hence, at freeze-out their multiplicity should be smaller 
than that of the a-particle. However, a more quantita- 
tive analysis is certainly needed. Within the equilibrium 
scenario the most stable nucleus might be Fe. However, 
little is known about the properties of Fe at finite temper- 
ature. At low densities one might expect a linear depen- 
dence of the binding energy (perturbative theory), how- 
ever the calculation for the a-particle clearly shows that 
this might not be valid for a stronger bound system. On 
the other hand during (central) heavy ion collision the 
dynamical generation of heavy nuclei needs time due to 
many (binary) collisions in the system and such heavy el- 
ements might not recombine at all in such a evaporation 
scenario |12| . 

A microscopic treatment of more complex (light) nu- 
clei could be achieved, e.g., along the lines of |37]. So far, 
larger and heavier clusters are " hidden" in the large num- 
ber of nucleons present. While those are rather easy to 
incorporate in the ideal gas picture, a calculation includ- 
ing medium dependence such as self energy corrections 
and Pauli blocking is more elaborate. 



5 Conclusion 

We have shown that a systematic microscopic calculation 
provides strong changes in the equilibrium composition of 
clusters in nuclear matter. The changes are strong enough 
to invert the ordering of multiplicities at freeze-out com- 
pared to the ideal case. Therefore an explanation of exper- 
imental results in terms of a microscopic picture with re- 
alistic nucleon nucleon forces evaluated in an equilibrium 
scenario might be possible. Note that a detailed compar- 
ison of this approach to the experimental data as, e.g., 
given by the INDRA collaboration for the SMM |9TTlT)] . 
needs much further investigation and has to be postponed 
to a future communication. 
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